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ABSTRACT 

With planets orbiting stars, a planetary mass function should not be seen as a low-mass ex- 
tension of the stellar mass function, but a proper formalism needs to take care of the fact that 
the statistical properties of planet populations are linked to the properties of their respective 
host stars. This can be accounted for by describing planet populations by means of a differen- 
tial planetary mass-radius-orbit function, which together with the fraction of stars with given 
properties that are orbited by planets and the stellar mass function allows to derive all statistics 
for any considered sample. These fundamental functions provide a framework for comparing 
statistics that result from different observing techniques and campaigns which all have their 
very specific selection procedures and detection efficiencies. Moreover, recent results both 
from gravitational microlensing campaigns and radial-velocity surveys of stars indicate that 
planets tend to cluster in systems rather than being the lonely child of their respective parent 
star. While planetary multiplicity in an observed system becomes obvious with the detection 
of several planets, its quantitative assessment however comes with the challenge to exclude 
the presence of further planets. Current exoplanet samples begin to give us first hints at the 
population statistics, whereas pictures of planet parameter space in its full complexity call 
for samples that are 2—4 orders of magnitude larger. In order to derive meaningful statistics 
however, planet detection campaigns need to be designed in such a way that well-defined 
fully-deterministic target selection, monitoring, and detection criteria are applied. The prob- 
abilistic nature of gravitational microlensing makes this technique an illustrative example of 
all the encountered challenges and uncertainties. 

Key words: planetary systems — gravitational lensing. 



1 INTRODUCTION 

More than 450 planets orbiting stars other than the Sun have been 
detected to date by means of four different techniques: Doppler- 
wobble stellar radial-velocity measurements, planetary transits, 
gravitational microlensing, and the direct detection of emitted or 
reflected light. Observing campaigns now need to evolve from 
the pure detection of planets to studies that allow to infer the 
statistical properties of the underlying populations that are being 
probed. In order to achieve such a goal, deterministic procedures 
for the selection of targets and the identification of planetary sig- 
nals are required jDomimk et al.ll200ll2008l ; iTsapras et ai]|2009t 
lO'Toole et al.ll2009HDominik et alfeoioh . 

While many studies on planet populations based on data from 

radial - velocity surveys have been carried out (e .g. Marcy et al. 

Udrv & Santos! l2007h iMarcv et al.l 120081 ; iMavor & Udry 



2005; 



2008: IO'Tooleetal.l2 009). there has been a long silence on extract 



ing planetary abundances from gravitational microlensing cam- 
paigns since the twin papers on the first five years of the PLANET 



campaigrQ jAlbrow et all 1200 it Gaudi et alj|2002h . and those on 
the OGLE-I llTsapras et alj|2003h and the first year of the OGLE- 
III surve )Q (Snodg rasset"aU l2004). Only very recently, there ap - 
pears to be a sudden inflation iSumi et alj2010tlGould et al.l 2010). 

The various techniques currently used for studying planet pop- 
ulations (and any future ones as well) have their very own preferred 
regions of planet parameter space that they are sensitive to, while 
being blind to others. Consequently, these do not directly probe the 
planetary mass function, but some bits and pieces for which the 
detection efficiency of the campaign and the selection biases need 
to be determined carefully. Any quoted planet abundance needs to 
come with a complete description which region of planet parameter 
space it refers to, and how it has been averaged. In particular, one 
can anticipate substantial differences on whether one talks about 
Solar-type stars or M dwarfs, hot or cool planets, bulge or disk stars 
(with their different metallicities). 

In order to get around such difficulties, and ease the compar- 
ison between findings that arise from different campaigns and/or 
techniques, a general framework of differential planetary mass 
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functions with respect to their fundamental parameters is suggested 
in this paper. 

Rather than planets being distributed just randomly amongst 
stars, it seems that they tend to cluster in planetary systems. Such 
a conjecture is underpinned by recent observational evidence on 
both outer gas-giant planets dGaudi et al . 2008; Marois et alj2008h 
as well as Su per-Earths and planets with Neptune-class m asses in 
closer orbits dMavor & Udrvll2008l : iLo Curto et alj|20ld) . There- 
fore, the detection of planets in an experiment does not correspond 
to independent draws from its population, but the probability of 
their detection around a star that is known to host planets is larger 
than that of finding it around a randomly chosen star. 

While Sect.|2]presents a theoretical framework for describing 
planet populations in view of clustering in planetary systems and 
specific regions of interest or sensitivity, Sect. [3] provides rough 
estimates for the size of planet samples required to assess the fun- 
damental functions that decribe these populations. Sect. [4] is de- 
voted to planetary multiplicity, while Sect.[5]looks into planet abun- 
dance estimates arising from gravitational microlensing observa- 
tions, their uncertainties, and the involved challenges. Sect. [6] fi- 
nally concludes the paper with a short summary and outlook. 



2 FUNDAMENTAL FUNCTIONS DESCRIBING 
PLANETARY SYSTEMS 

Only celestial bodies that are in orbit around a around a star or 
stellar remnant are 'planets' This means that planets cannot be 
seen in isolation from these, and consequently planets are not well- 
descr i bed by just extending th e stellar mass function |Salpeter 
1 19551 : Scalo 19861; I Kroupall 2002) to lower masses, but a mass func- 
tion decribing planets needs to link to their host stars (or remnants). 

Let us consider explicitly the dependence of planetary abun- 
dance on stellar mass M*, metallicity Z, age t, and spin 
rate SI and therefore define a differential stellar mass function 
£(M*, Z, r, f2)Q so that for a population with density functions 
pz(Z), Pt(t), and pn(Sl) for metallicity, age, or spin rate, respec- 
tively, one obtains a mass function 

S(Af*)= / ^M i ,,Z,T,Sl)p Z (Z)dZp T (r)dTpn(Sl)dSl, (1) 



(2) 



x d[lg(r p /r©)]d[lg(a/a®)]de, 



(3) 



where the number density of stars becomes 



S(M*)d[lg(M*/M )], 



where Mq denotes the mass of the Sun. 

For the stars that host planets, the properties of planets can 
then be described by a differential planetary mass-radius-orbit 
function (p(m p , r p , a, e; M*, Z, r, SI), where m p , r p , a, and e de- 
note the mass, radius, orbital semi-major axis, or orbital eccen- 
tricity of the planet, respectively, and further parameters might be 
added. This implies a mass function for planetary systems around 
stars with (M*, Z, t, SI) given by 



$(m p ;Af*,Z, r,fi) 



(p(m p , r p , a, e; M*, Z, r, SI) x 



3 following the "Position Statement on the Definition of a 'Planet'" (in the 
revision dated 28 February 2003) by the IAU Working Group on Extrasolar 
Planets (WGESP) 

4 For the vast majority of stars, the spi n rate O essentially becomes a func- 
tion of stellar mass M* and age r (e.g. ICollier Cameron & Li ; 1994; Barnes 
2007), so that this parameter can be neglected. 



with ?-0 being the Earth's radius and a® = 1 au, so that the average 
number of planets in such systems reads 

n p {M„,Z, t, n) = J $(m p ; M 4 , Z, r, SI) d[lg(m p /M ffi )] , (4) 

where M® is the mass of the Earth. 

With / P (M*, Z, t, SI) denoting the fraction of stars that host 
planets, the number density of planets for a stellar population be- 
comes 



N P = J f P (M* , Z, r, Q) e(M* , Z, r, SI) n p (M* , Z, r, SI) x 

x p z {Z) dZp r {r) drpn(Sl) dSl d[lg(M/M*)] . (5) 

Moreover, one finds a population-integrated planetary mass- 
radius-orbit function 



ip(m p ,r p ,a,e) = J / P (M*, Z,r,Sl) £(M*, Z,r,Sl) x 
x ip(m p ,r p , a,s; M*, Z,t,SI) x 

x Pz (Z)dZp r {T)dTpn(Sl)dSld[lg(M*/M e )], (6) 
and a corresponding planetary mass function results as 



*(m p ) = y ^(mp,r p ,a,e)d[lg(rp/r-®)]d[lg(a/ae)]de, (7) 
so that one finds the number density of planets again as 

N p = [ *(m p )d[lg(m p /M ffi )]. (8) 



Provided that experiments in the hunt for extra-solar planets 
follow deterministic criteria, a mass function can be extracted that 
refers to the selected host stars and planetary orbits that the ap- 
plied technique is sensitive to, i.e. averages are taken over the stel- 
lar population and the orbital parameters. However, in order to an- 
swer fundamental questions such as 'How frequent are planets of a 
given mass range in the Solar neighbourhood?', 'What fraction of 
stars in the Milky Way do have planetary systems?', or 'How many 
planets that could host life are there in the Universe?', one needs 
to trace back the description of planetary systems to more funda- 
mental functions such as the differential mass-radius-orbit function 
<p(m p , r p ,a, e; M*, Z, t, SI), the fraction of stars with planetary 
systems / p (Af*, Z, t, SI), and the differential stellar mass function 
t{M*,Z,r,Sl). 



3 MEASURING PLANETARY MASS FUNCTIONS 

In order to obtain an estimate on how well we can measure a plan- 
etary mass function, let us consider dividing the parameter space 
into multi-dimensional bins. If rather than aiming for a precision 
measurement, one sets the goal at an 'astronomical' accuracy of 
50 per cent, the assumption of Poisson statistics yields the re- 
quirement of each bin to contain at least 4 planets. Let p denote 
the number of considered parameters, and b the number of con- 
sidered parameter ranges, the minimal number of planets needed 
to provide the desired result is N p = 4 6 p , which would cor- 
respond to letting the choice of parameter ranges follow the ob- 
served distribution of detected planets, in such a way that each bin 
contains exactly the mimimum of 4 planets. With k denoting the 
desired accuracy, one finds more generally iV p = k~ 1/2 b p . Ta- 
ble Q] shows the requirements for some selected cases with rela- 
tive accuracies of 50 per cent or 20 per cent, 2-, 4- or 6-parameter 
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functions, and a various number of bins ranging from 2 to 10. 
Roughly, one gets an idea of the distribution of the planet abun- 
dances with 6^3, but one can realistically only start talking about 
a "planetary mass function" for & ^ 5. While a planetary mass- 
radius-separation function ipm p ,r p ,a(m p , r p , a; M*, Z, r) depend- 
ing on the stellar mass, metallicity, and age involves 6 param- 
eters, less detailed 4-parameter functions are e.g. the planetary 
mass-separation tp m p ,a(m p , a; M*, Z) or mass radius function 
fm p .j-p (m p , r p ; M* , Z) depending on stellar mass and metallicity, 
or a planetary mass-radius-separation function depending on stel- 
lar mass only, and 2-parameter functions would e.g. be the plane- 
tary mass function <fim p {m p ; M*) depending on stellar mass only, 
or the planetary mass-separation function tp m a (m p , a) irrespec- 
tive of the stellar properties. We now have a total sample of about 
450 planets orbiting stars other than the Sun, where it took about 
10 years to detect the first 150, then about 3 years to detect the 
next 150, and then just about 1 year to detect the equal number 
of 150. TableQ]shows how long campaigns with a constant detec- 
tion rate of 150 planets per year would have to last in order to obtain 
the respective functions with desired accuracies. 

Right now, the collected data allow to measure 1 -parameter 
functions, find the basic structure structure (b 10) of 2-parameter 
functions, see basic trends (b ^ 3) in 4-parameter functions, and 
some hint on the dependency of the planet abundance on further 
parameters. With 150 planets per year, or more realistically, a fair 
factor of this rate, rough ideas (b ^ 5) of 4-parameter planetary 
mass functions (b ^ 5) and an indication of trends (b ^ 3) for 
6-parameter planetary mass functions are obtainable within fore- 
seeable time frames, but the numbers call for more aggressive 
searches. 



4 PLANETARY MULTIPLICITY 

While stars with and without planets have been distinguished by 
referring to the fraction / P (M*, Z, t, fl) of stars that host planets 
and defining the differential planetary mass-radius-orbit function 
<fi(m p , r p , a, e; M*, Z, t, fi) to relate to these only, a further statis- 
tic is the distribution of the number of planets amongst all plane- 
tary systems. With multiplicity indices (k that denote the fraction 
of planetary systems containing k planets, where 

oo 

£Ck = l, (9) 

fc=i 

the planetary mass-radius-orbit function can be decomposed as 

p(m p , r p , a, e; M*,Z, r, Q.) 

OO 

= ^fcCfc^fc(m p ,r p ,a,e;M*,Z,r,n) , (10) 
fc=i 

where 

J <p k (m p ,r p ,a,E]Mi<,Z,T,£l) x 

x d[lg(m P /M®)] dpg(r P /r©)] dpg(o/o e )] Ae = k . (11) 

In general, all ipk(m p ,r p , a, e; M*, Z, t, fi) may be different. 
Together with the multiplicity indices £fe, one would be left with an 
infinite number of parameters. This however can be meaningfully 
avoided by adopting a functional dependence of and ipk on k 
that is described by a small finite number of parameters. 

In particular, one might want to distinguish stars with a single 
planets to multiple-planet systems, described by (with £ mu it = 
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N p 
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0.5 
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4 x 5 6 = 
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420 years 
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4 x 10 6 = 


= 4, 000, 000 


27,000 years 
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25 x 2 2 = 
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8 months 
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2 


3 


25 x 3 2 = 
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1.5 years 
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2 
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25 x 5 2 = 
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4 years 
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10 


25 x 10 2 
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17 years 
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25 x 2 4 = 
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25 x 3 4 = 
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13.5 years 


0.2 
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25 x 3 6 = 


= 18,225 


120 years 
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25 x 5 6 = 
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10 
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= 25, 000, 000 


170,000 years 



Table 1. Minimal number of planets N p required to sample a descriptive 
statistic with p parameters with b bins to a relative accuracy re, and time 
T150 required to acquire such a sample for a planet detection rate of 150 
per year. 



1 - CiX 0i(m py r p ,a,e; M+,Z,t, Q), and 
<£mult(m p , r p , a, e; M*, Z, r, Q) 

OO 

= ^^kC,k<Pk{m p ,r p ,a,E;M*,Z, t,Q) 

k = 2 

= ip(m p ,r p ,a,£\ M*, Z,t,Q) — 

~ Ci 0i ( m p> r pi a i e i M+, Z, r, Q) . (12) 

In fact, IWright et alj j2009h have argued that there is evidence for 
if i being different from mu it. 

The assessment of planetary multiplicity however poses a 
huge challenge for properly interpreting the observational data, 
given that our knowledge of the absence of further planets in ob- 
served systems is quite limited. If Hot Jupiters are considered 
lonely, whereas Neptune-mass planet s are frequently found i n mul- 
tiple systems jMavor & Udrvl 120081 ; ILo Curto et at] |201Ch . how 
much does this have to be attributed to the fact that observational 
techniques that report Hot Jupiters are insensitive to less massive 
planets, whereas if the sensitivity extends down to lower masses, 
other such planets are spotted rather easily? It is intriguing to see 
that observations of transit timing variations led to the suggestion 
of the presence of a 15 E arth-mass planet in the WASP-3 system 
dMacie iewski et al. 20 H^) th at was already known to host a Hot 
Jupiter JPollacco et alj|2008l) . Planets reported by microlensing in 
particular cannot be claimed to be the only ones in the system, they 
were just th e only ones t hat rev ealed their presence during a tran- 
sient event. iKubas et alj d2008t) explicitly found that the acquired 
data do not exclude the presence of gas-giant planets at any separa- 
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tion orbiting the lens star that caused event OGL E-2005-BLG-3 90, 
which is known to host a cool Super-Earth teeaulieu et alj2006t) . It 
is easier to detect a planet than being able to claim that there are no 
other planets orbiting the same star, and if one aims for quantifying 
multiplicity, this needs to be addressed. 

Planetary multiplicity however becomes an obvious phe- 
nomenon with the detection of respective systems, suc h as the pair 
of ga s-giant planets orbiting OGLE-2006-BLG-109L l lGaudi et al.l 
2008), which resemble a half-scale version of the Jupiter- Saturn 
part of the Solar system. Interestingly, the planets found to or- 
bit HR 8799 look like the complementary double-scale version 
dMarois et al]|2008h . It is particularly striking that very early op- 
portunities to detect such systems by gravitational microlensing or 
direct imaging, respectively, were successful, while one needs to 
keep in mind that planets with an orbital period similar to Saturn 
cannot be detected from radial- velocity surveys so far (given a 10- 
15 year history of respective campaigns), and observing planetary 
transits is further disfavoured by the small transit probability. How- 
ever, for Super-Earths and planets with Neptune-class masses in 
closer orbits, rad ial- velocity surveys find a very high level of mul- 
tiplicity as well dMavor & Udry||2008l ; lLo"Curto et alj|2010t) . 

The detection of the pair of Jupiter- and Satu rn-like planets 
orbiting OGLE-2006-BLG-109L jGaudi et alj2008l) is often hailed 
because of the striking similarity with the Solar system, albeit that 
there is basically nothing that can be said about potential inner 
rocky planets other than that such cannot be excluded. There is 
however another important result arising from this discovery: outer 
gas-giant planets are not of the lonesome type. How does one arrive 
at such a conclusion? Regardless of the large detection efficiency 
for such planets in events with a peak magnification as large as that 
of OGLE-2006-BLG-109 (A ~ 290), the planetary abundance is 
moderate or small. If we consider an abundance of 5 per cent, the 
probability for a double catch would be just 0.25 per cent if the 
planets were drawn independently from the population. This would 
mean an expected detection of ~ 1/30 systems am ongst the 13 
events comprising the systematic sample reported by Gould et ajj 
d2010h . so that we would have been very lucky to find the detected 
pair. Therefore, it appears the more likely assumption that the two 
detections were not the result of independent draws, but instead the 
probability for a planet to orbit a star is larger if one considers a star 
that is known to host planets as compared to an arbitrarily chosen 
star that might host planets or not. This however means that it is 
not appropriate to consider a planetary mass function with planets 
randomly drawn from it, but instead one needs to distinguish be- 
tween stars with or without planets, as the formalism suggested in 
the previous section does. These arguments however get weaker if 
the planetary abundance was as large as 20 per cent, because this 
would mean a probability of 4 per cent for a pair, or 1/2 expected 
to be detected for 13 events as compared to the one found. 



5 PLANET ABUNDANCE ESTIMATES FROM 
MICROLENSING OBSERVATIONS 

The statistical analysis of microlensing events in order to derive 
planet abundance estimates provides an illustrative example of the 
challenges one is facing. If a planet orbits the lens star, a detectable 
signal will only arise with a finite probability. This finite detection 
efficiency for planets of given mass and orbital separation from 
their host star is of relevance not only for assessing abundances 
by means of detections but also for drawing conclusions from the 
absence of planetary signals. Moreover, the host stars of planets 



detected by gravitational microlensing arise stochastically from the 
underlying population of stars that intervene the observed targets, 
with current experiments most of the masses of the lens stars are 
only known up to a broad probability distribution, although the 
mass of the lens star is frequently known for events i n which plan- 
etary signals have been detected dGould et alj|2010l) . The lack of 
information about the planet's host star is troublesome, since it is 
important to distinguish planet population statistics in the range of 
stellar masses between 0.1 and 0.8 Mq, covered by microlens- 
ing, given that current planet-formation models predict substan- 
tial differences, i n particular for the abundance of gas-giant planets 
dlda & Linl2005r) . 

For the rather small planet samples acquired so far, let us 
however neglect this issue for the time being, and just compare 
planet abundance estimates that refer to the sample of probed lens 
stars. Recently, there have been some discussions about planetary 
mass functions that can be e xtracted from micr olensing observa- 
tions. While the discussion bv lSumi et all fcoiO) is not based on a 
well-defined criterion for selecting the considered 10 planet detec- 
tions f rom the so far p ublished 24 candidates towards the Galactic 
bulge toominikl feoiO). and moreover no relation has been given 
between th ese 'detections' and the efficiency of the full observing 
campaigns, iGould et alj J20101) in contrast adopted selection cri- 
teria that lead to a well-defined event sample, and evaluated the 
detection efficiencies properly. However, they refer to a planetary 
mass function described by means of the planet-to-star mass ratio, 
whose value is highly questionable, given that rather obviously one 
does not expect the same number of half-massive planets to form 
around half-massive stars. In particular, coagulation and accretion 
processes depend on the masses of the bodies involved and their 
spatial density, b ut not on the mass of the star. Nevertheless, the 
sample drawn bv lGouldetallfeOlOl) allows for an insightful fur- 
ther look. 

IGould etail j201Ch refer to 13 events with a peak magnifica- 
tion Ao > 200, densely monitored by MicroFUN_| (and other cam- 
paigns) from 2005 to 2008. Amongst those events, 2 provided a sig- 
nal that indicates the presence of a massive gas-giant planet above 
150 M m (0.5 Mjup), namely OGLE-2006-BLG-109 and MOA- 
2007-BLG-400. For such planets, the detection efficiency for or- 
bital separations that correspond to the 'lensing zone can broadly 
be as sumed to be of the order of 100 per cent dGriest & SafizadeH 
1998). One would therefore estimate the abundance of such planets 
to be about 15 per cent. 

Rather than just focussing on event s with large peak 
magnifications, the PLANET collaboration jAlbrow et al.l 1 199*81 : 
Dominik et al. 2002) has acquired data on a much larger sample 
of about 50 events per year with Ao > 2 from 2002 to 2007, with 
sampling intervals of around 2 hrs or better, where the average de- 
tection efficiency for Jupiter-mass planets i n the 'lensing zone' fo r 
such a sample is about 15 to 20 per cent dGould & Loeb|[l992l) . 
Only o ne respective planet has been reporte d: OGLE-2005-BLG- 
071Lb dUdalski et al Jl 20051 : IPong et aljboO^) . as compared to ex- 
pected 45-60 if those reside around each of the lens stars. This 
gives a rough abundance estimate of 1.5-2 per cent, which looks 
substantially smaller than what one guesses from the densely mon- 
itored events with Aq > 200. The PLANET team earlier claimed 

5 http://www.astronomy.oWo-state.edu/~microfun/ 

e Angular planet-star separations 0.618 £?e ^ ^ 1.618 where 
6 E = {[(4GA/)/c 2 ] (D' 1 - D~ 1 )} 1/2 , with D L and D s the distances 
of the lens and source star from the observer, respectively, typically a range 
of 1.5 to 4 au. 
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an upper abundance limit (at 95 per cent confidence) of 33 per cent 
on Jupiter-mass planets in the same orbital range based on the ab- 
sence of any detection amongst 42 events well -covered from 1995 
to 1999 jAlbrow et alJl200ll ; lGaudi et alj|2002l) . 

The results based on the MicroFUN and PLANET data how- 
ever appear to be statistically compatible, and one finds that 
the small-number statistics imply large uncertainties. In fact, an 
Agresti-Coull confidence interval for the plan etary abundance 
based on the underlying binomial distribution dAgresti & Coulll 
Il998h at 95 per cent probability extends from 3 per cent to 43 per 
cent for the MicroFUN result, and from less than 0.01 per cent to 
10 or 13 per cent for the PLANET result. This also gives some in- 
dication o f the unqual ified uncertainties of the statistical results 
quoted bv lGould et ail feOld) . 

Apart from the binomial statistics, there are some system- 
atic uncertainties. One might in fact wonder whether there are 
further planets just waiting to be detected in the PLANET data 
that have not been spotted yet due to absence of a comprehen- 
sive systematic analysis. This may not be too unlikely, given that 
e.g. the event OGLE-2008-BLG-513 was initially considered to 
be due to a stellar binary, before a planetary model had been sug- 
gested jGould et al1l2010h . On the other hand. lCohen et alj 12010) 
have pointed to a puzzle regarding the properties of the high- 
magnification events that casts doubt on whether we really under- 
stand the mechanism responsible for producing these. Namely, a 
highly statistically significant correlation has been found between 
the event peak magnification and the metallicity of the observed 
source star (i.e. not the planet's host star). This sample bias is so far 
not understood. 

If the planetary abundance turns out to be small, its determina- 
tion becomes more difficult. If one focuses on high-magnification 
peaks, for an abundance 15 per cent, 5 per cent, or 2 per cent, the 
monitoring of 90, 320, or 800 microlensing events, respectively, 
would be required in order to make the half-width of a symmetric 
95 per cent confidence interval match half the abundance. With an 
average detection efficiency of 15-20 per cent for less favourable, 
but useful {Aq ^ 2), hourly-sampled events one would require 5- 
6 times as many events, but given that the number of events with 
^4.0 > 200 is about 100 times smaller, such a strategy looks more 
feasible, in particular since with the current detection rate of the 
microlensing surveys, the monitoring of ~ 200 suitable events per 
year is possible. 

While the detection efficiency for Jupiter-mass planets is a 
rather robust number, the sensitivity of microlensing campaigns to 
planets between 1 and 10 Af® ("Super-Earths") is a strong function 
of planet mass and orbital separation. Therefore, the interpretation 
becomes substantially dependent on the choice of the considered 
region of planet parameter parameter space and the adopted averag- 
ing. The sparcity of data makes a meaningful assessme nt quite diffi- 
cult. T he least massive planet in the sample adopted bv lGould et al.l 
and the only one below 50 Af®, was found to have a mass 
of 13 A/® with a substantial uncertainty, so that it may or may 
not fall into the Super-Earth mass regime. Moreover, detection ef- 
ficiencies in this region vary substantially amongst the 13 events 
that comprise the sample with a prominent peak close to the an- 
gular separation of the planet being equal to the angular Einstein 
radius of its host star. The situation is slightly better for the results 
of the PLANET campaign: OGLE-2005-BLG-390Lb has been es- 
timated to h ave a mass between 3 and 10 Af® w ith a probability of 
68 per cent dBeaulieu et al ,l2006l : lDomini32 006), and one is within 
the right order of magnitude by assuming a detect ion efficiency of 
about 1-2 per cent for 'lensing zone' planets (c.f. lBennett & Rhiel 



1 19961 : 1 Kuba s et al ] |2008h for an average over 300 observed events. 
One detection then leads to an maximum-likelihood point estimate 
for the respective abundance of 17-33 per cent, about 10-20 times 
larger than that obtained for planets above 0.5 Afj up from the same 
set of observations. However, none of the derived estimates should 
be considered to be correct to within less than a factor 3^1, but 
on the other hand, they are not substantially worse either. While 
microlensing observations show that it is implausible that low- 
mass planets are rare, drawing firm conclusi ons is prevented by 
the current low-number statistics. In particular. iGould et alj J201CI) 
found that their sample size is insufficient for reliably determining 
a power-law index in the mass function with respect to the planet- 
to-star mass ratio. However, they estimate the local planet number 
density in orbital separation d = 9 p /9e and mass ratio q at d — 1 
and q = 5 x 10~ 4 to be 0.36 per decade in each of the quantities. 
With power laws considered to range from oc q~ ' 2 to oc q~~ o e , and 
assuming a 'typical' stellar mass of 0.3 Af®, one finds for the mid- 
dle of the decade from 1 Af® to 10 Af®, roughly at m p = 3 Af®, 
local number densities of 0.6 or 1.9, respectively. For comparison, 
let us account for the fact that the lensing zone covers 0.42 decades, 
and therefore multiply the derived abundance of 17-33 per cent by 
2.4, which resu lts in values 0.4-0 .8, not much different from the 
result found bv lGould et ail ( l2010h , where one also needs to con- 
sider that the local density is not equal to the average density in 
the considered surrounding region. Applying the sa me procedure 
for tra nsferring the local number density according to Go uld et all 
(2010) to a planet mass m p = 1 Mj up yields values of 0.25 or 0.12, 
for the two different power-law indices, respectively. The range of 
power laws gives abundance density ratios between m p = 3 Af® 
and m p = 1 Mj up in the range from 2.5 to 16 (as compared to 
10-20 guessed from PLANET observations). 



6 SUMMARY AND OUTLOOK 

A planetary mass function is not the extension of the stellar mass 
function to lower masses, given that planets orbit stars. Instead, 
planetary mass functions more or less strongly depend on the char- 
acteristic proper properties of the respective host stars, such as the 
stellar mass Af*, metallicity Z, age r, and spin rate f2 as well. As 
long as the planet mass m p , planet radius r p , orbital semi-major 
axis a, and orbital eccentricity s are considered as descriptive pa- 
rameters, all population statistics can be derived from three fun- 
damental functions, namely the differential mass-radius-orbit func- 
tion </s(m p , r p ,a,e; Af*, Z, t, fl), the fraction of stars with plan- 
etary systems / p (Af*, Z, t, fl), and the differential stellar mass 
function £(Af*, Z, r, f2). In principle, tne fragmentation of the 
planetary system into k planetary bodies gives us an infinite number 
of multiplicity indices £fc that correspond to the fraction of plane- 
tary systems with exactly k planets, as well as respective specific 
mass-radius-orbit functions tf>k(m p , r p , a, e; Af*, Z, r, fi). Adopt- 
ing a functional dependence on k however allows for a descrip- 
tion with a finite number of parameters. A first-order step would 
be to distinguish single-planet and multiple-planet systems. How- 
ever, the determination of their respective fraction for all stars that 
host planets requires a proper assessment of the constraints on the 
presence of further planets in the studied systems. 

It cannot be stressed enough that the formation and evolution 
of planets is a crucial step towards the development of life, but it 
will not be understood by focusing the interest on habitable plan- 
ets, rather than embracing planet populations in their amazing di- 
versity. Moreover, rather than just optimizing planet searches for 
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a large detection rate, the critical design feature is to follow well- 
defined monitoring and detection criteria that allow to carry out 
simulations, so that meaningful statistics can be derived. 

With the probabilistic nature of the alignment of stars suit- 
able to provide gravitational microlensing events and similarly for 
a planet to reveal its presence around the foreground (generally 
unobserved) lens star, studying planet populations by microlens- 
ing prominently comes will all difficulties and challenges that one 
might encounter. All conclusions that can be drawn to date not 
only suffer from small-number statistics, but moreover from diffi- 
culties in the assessment of the detection efficiency and the related 
need to refer to strictly deterministic procedures for the monitoring 
strateg y (e.g. lDominik et al.ll2010l : [Gould et alj|2010h . lGould etaU 
(2010) have recently presented a first statistically meaningful, but 
very small, sample of 13 events comprising only those events for 
which source and lens stars are so closely aligned to yield a peak 
magnification Aq ^ 200. While some fundamental statistics can 
be derived to within a factor 3-4, current data do not allow to 
do substantially better. The fact that a maximum-likelihood point 
estimate for the abundance of planets above 0.5 Mj up from this 
sample comes out as 10 times larger than what one would guess 
from 6 years of PLANET observations that include events with 
smaller peak magnifications as well, while both values are statisti- 
cally compatible, demonstrates the current uncertainties. Moreover, 
one might wonder whether any systematics that are not fully under- 
stood affect this outcome. The rel atively large numbe r of gas-giant 
planets in the sample adopted by I Gould et alj d20 lOh as compared 
to the PLANET observations with an about 3-5 times larger to- 
tal detection efficiency is somewhat surprising, in particular given 
that such p lanets are readily detected at high efficiency already for 
A 10 dGriest&Safizadeh| [l998). In fact, for measuring plan- 
etary abundances, focussing on the few events with very large Ao 
does not appear to be a promising strategy, due to the rarity of such 
events. 

The prospects for obtaining a measurement of the planetary 
abundance with a desired accuracy strongly depend on the abun- 
dance itself, given that a small abundance will imply a small num- 
ber of detections. This is serious limiting factor for any planet de- 
tection campaign, and the measurement of the abundance of 'sec- 
ond Earths' by NASA's Kepler mission is not immune to the prob- 
lem of small-number statistics either. Rather than trying to esti- 
mate an abundance from a small region near the sensitivity limit, a 
more robust estimate for the abundance of habitable planets would 
arise from an interpolation between hotter and cooler planets, mak- 
ing use of a larger number of detected objects, and assuming that 
planet formation will not be radically different just for the habit- 
able zone. The improvement of the statistics by a more powerful 
mission such as PLATO (PLAnetary Transits and Oscillations!] 

i — | -n 

(Catala et al. 2009) appears to be much desired. 

The dependency of the planet population statistics on the prop- 
erties of the host star implies that a quite substantial amount of data 
will need to be collected for properly measuring the planetary mass 
function and even more for determining the full mass-radius-orbit 
function. With the 450 reported planets so far, we can now assess 
1- and 2-parameter functions (if and only if we understand any de- 
tection bias), but it requires a much larger detection rate than the 
recent 150 planets per year in order to understand the distribution 
of planets in the Universe. From 1995 to 2009, the respective time 
interval of acquiring 150 planets has been cut by a factor three twice 
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(10 years, 3 years, and then 1 year), so that we see a substantial ac- 
celeration, and the Kepler mission is already contributing to boost- 
ing the planet detection rate further. Nevertheless, drawing pictures 
of planet parameter space in its full complexity calls for samples 
that are 2-4 orders of magnitude larger than those we have now. 
Looking back at the history of exoplanet detections however tells 
us something else: we gained a lot of insight from a few individ- 
ual detections that came as a surprise and challenged the prevailing 
understanding. Shouldn't we expect to be surprised again when em- 
barking on exploring further uncharted territory? Our understand- 
ing of planetary formation and evolution is not probed uniformly 
by a planetary mass-radius-orbit function, but certain regions of 
planet parameter space may prove more critical in the power to 
discriminate between alternative theories or to measure crucial pa- 
rameters. Therefore, efficiency could be gained from observational 
campaigns delivering specific characteristic statistics that can be 
robustly determined from rather small samples. 
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